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We explore from a theoretical perspective angular distributions of electrons emitted from a Na§ cluster after excitation 
by a short laser pulse. The tool of the study is time-dependent density-functional theory (TDDFT) at the level of the 
local-density approximation (LDA) augmented by a self-interaction correction (SIC) to put emission properties in 
order. We consider free Na§ and Na§ deposited on the surfaces MgO(OOl) or Ar(OOl). For the case of free Nas, we 
distinguish between a hypothetical situation of known cluster orientation and a more realistic ensemble of orientations. 
We also consider the angular distributions for emission from separate single-electron levels. 



1 Introduction Optical methods have provided the 
key analyzing tools in cluster physics over decades. In 
the first stage, optical absorption measurements allowed 
one to collect rich information on structure and dynamics 
of these small, nano-sized particles; for an overview see, 
e.g., IBS!. More information can be gathered when one 
additionally measures the features of reaction products in 
more detail. A most important channel in this context is 
electron emission and thus there meanwhile exist many 
investigations that analyze the properties of the electrons 
emitted after irradiation by a short laser pulse. The first 
step in that direction is photoelectron spectroscopy (PES) 
where the distribution of the kinetic energy of the emitted 
electrons is recorded. This method has been applied in 
cluster physics since long, see e.g. |S]|6]. Stepping further 
in refinement, one can determine the angular distribution 
of the outgoing electrons which, in fact, is mostly done 
simultaneously together with PES 171 4121 . Further chal- 
lenging aspects come into play when considering clusters 
in contact with substrates. That combination is often moti- 
vated (if not dictated) by easier experimental handling. It is 
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of great importance for possible practical applications, and 
the effects at interfaces are also an interesting problem for 
basic research, for up to date collections see e.g. II 1311 141 . 
It is obvious that the deposition on a surface modifies the 
emission properties, in particular angular distributions, 
thus calling for dedicated theoretical studies. This paper 
is devoted to a theoretical exploration of angular distribu- 
tions of laser excited metal clusters, free and deposited on 
insulating surfaces, MgO(OOl) and Ar(OOl). 

There is a broad choice of approaches for the descrip- 
tion of clusters on surfaces, from fully detailed quantum 
mechanical treatments of all constituents lfl5ll to a robust 
jellium treatment of cluster and interface |[T6l . We are go- 
ing here for an intermediate strategy, detailed at the side of 
the highly reactive cluster and less so for the inert substrate. 
As a basis of the description, we employ density-functional 
theory ifTTl . which provides a most efficient theoretical de- 
scription for the electronic structure and dynamics of clus- 
ters lfl8l . We simulate the detailed dynamics of laser exci- 
tation and subsequent electron emission by time-dependent 
density-functional theory (TDDFT) solved on a grid in co- 
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ordinate space. The substrates are inert and will be de- 
scribed classically in terms of polarizable atoms I119II20I . 
The computation of angular distributions requires rather 
large numerical boxes. Therefore, previous explorations 
dealt with approximations, a semi-classical approach [21 1 
which is confined to very intense laser pulse or a quantum- 
mechanical, but two-dimensional axial, treatment 12211231 
for lower intensities and short pulses. Here we are going 
for a fully three dimensional TDDFT analysis because both 
simplifications are not applicable anymore for the studies 
intended here and because the progress of computational 
resources meanwhile allows a TDDFT analysis in full three 
spatial dimensions. We will take Nag as a test case and 
compare free Nag with Nag on MgO(OOl) as well as Nag 
on Ar(OOl). 

The study concentrates on direct electron emission, 
i.e. those electrons which are leaving the cluster with or 
directly after laser impact. In practice, there is a com- 
petition between direct and thermal emission. The pre- 
ferred exit channel depends on the length of the laser pulse. 
For pulses longer than the collisional relaxation time in- 
duced by two-electron processes, thermal emission domi- 
nates, while shorter pulses change the weight towards di- 
rect emission 124-271. We will use here pulses with full 
width at half maximum of 20 fs, which surely are at the 
side of dominant direct emission. 

There is also a subtle point about the orientation of 
the cluster relative to the laser polarization. Experiments 
with clusters in the gas phase deal, in fact, with an equi- 
distributed ensemble of orientations, while deposited clus- 
ters have a well defined orientation due to the known sur- 
face structure. We will briefly address this question, con- 
sidering both free clusters with known orientation and av- 
eraging over cluster orientations. We will then return to ori- 
ented clusters and investigate the changes caused by depo- 
sition on a surface. We will discuss how the specific angu- 
lar distributions for emission of each single-electron state 
separately combine to the total distributions and how this 
helps to understand the pattern observed in emission from 
deposited clusters. Furthermore, we will investigate and 
compare two different substrates, MgO(OOl) and Ar(001) 
surfaces. Both are insulators with a rather large band gap. 
The MgO(OOl) surface shows more corrugation and has 
a stronger interface attraction, while Ar(001) is softer and 
dominated by core repulsion. 

The paper is outlined as follows. In section [2] we 
briefly summarize the formal framework, numerical strate- 
gies, and the basic properties of the test cases. In section 
[3] we present and discuss results on free Nag and also ex- 
plore the double differential distributions. Section H] com- 
pares results for free Nag to those for Nag deposited on 
MgO(OOl) and Ar(001). Conclusions are given in section 

El 

2 Brief summary of the model 



2.1 The degrees of freedom The hierarchical 
quantum-mechanical-molecular-mechanical (QM/MM) 
model has been described in detail elsewhere (see J20 28 ] 
and l29l for a recent detailed review). Here we give a brief 
outline. The constituents of the system and their degrees 
of freedom are: ip n (r,t) , n = 1, N e \ for valence elec- 
trons of the Na cluster, R j(N a) , i (Na) = 1, ...,Ni for the 
positions of the Na + ions, R^c) , = 1, ...,M for the 
positions of the O cores, R,-(«) , v- v > = 1, M for the 
centers of the O valence clouds, R^r.*) , i^ k ' = 1, M for 
the positions of the Mg 2+ ions. The Na cluster is treated in 
standard fashion: The valence electrons quantum-mecha- 
nically and the ions classically M18II30I . The MgO substrate 
is composed of two species : Mg 2+ cations and 2 ~ an- 
ions. The cations are electrically inert and can be treated 
as charged point particles; they are labeled by v- k >. The 
anions are easily polarizable; this fact is accounted for by 
associating them with two constituents : A valence elec- 
tron distribution (labeled by i^">) and the complementing 
core (labeled by i^). All ions of the MgO substrate are 
described as classical degrees of freedom in terms of po- 
sitions Retype) . The difference R( c ) — R^ describes the 
electrical dipole moment of an 2 ~ anion. The Mg and 
O ions are arranged in crystalline order corresponding to 
bulk MgO. The dynamical degrees of freedom for Mg and 
O are taken into account in an active cell of the MgO(OOl) 
surface region underneath the Na cluster. The active cell 
is embedded in an "outer region" of MgO material where 
only ions are kept fixed, while oxygen dipoles remain 
fully dynamical. Beyond that region, only the Madelung 
potential is considered. The effect of the outer region on 
the active part is given by a time-independent shell-model 
potential taken over from OTI . Actually, the substrate con- 
sists of six layers each containing 784 Mg 2+ ions and 784 
2 ~ ions. The ions in the lowest layer are fixed to pre- 
vent them from relaxing and forming an artificial second 
surface. The active cell consists of three layers, each con- 
taining square arrangements of 242 Mg 2+ cations and 242 
O 2- anions. 

In the case of Ar, the modeling follows a similar, al- 
though simplified, track 11911321 . The Ar substrate is com- 
posed of only once species, neutral Ar atoms, each of 
which being characterized by its position and dipole mo- 
ment (exactly as O 2- anions). The substrate comprises 384 
atoms; this model has been shown to provide a fair approx- 
imation to bulk I33II34I . 

2.2 Energy and fields The total energy is decom- 
posed as E = £?Na + -Esurf + -^coupi where E-^ a describes 
an isolated Na cluster, .E Surf the MgO(OOl) or Ar (001) 
substrate, and E coup \ the coupling between the two subsys- 
tems. For En a , we take the standard functional of TDDFT 
at the level of the local density approximation (LDA), cou- 
pled with molecular dynamics (MD), as in previous studies 
of free clusters I18II30II including an average density self- 
interaction correction l35ll . The energy of the substrate and 
the coupling to the Na cluster consists of the long-range 
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Coulomb energy and some short-range repulsion which is 
modeled through effective local core-potentials B1I36I37I . 
To avoid the Coulomb singularity and to simulate the fi- 
nite extension of Ar, Mg 2+ and O 2 " ions, we associate 
a smooth charge distribution p(r) oc exp(— r 2 /<r 2 ) with 
each of these ionic centers. This yields a soft Coulomb po- 
tential to be used for all active particles. The model param- 
eters were calibrated to represent results of calculations 
where a fully quantum-mechanical description of the ac- 
tive MgO cell was employed in the case of MgO, for de- 
tails see EOl . The model parameters for the Ar substrate 
were chosen to reproduce model calculations of the NaAr 
dimer [38 1, in turn fitted to ab-initio data. 



Ai(r) defined as : 



The field equations obtained by variation of the above 
energy are augmented by the external laser field 



^las = -er-np ol i;osin 2 



t 



pulse 



7T I sin (wia S t) (1) 



which is activated only in the time interval < t < 
Tpuisc- The field strength E is related to the intensity / 

as l/(Wcm- 2 ) = 27.8 [fib/fVcm" 1 )] 2 . The total pulse 
length Tp U i sc corresponds to a full width at half (intensity) 
maximum as FWHMw T pu i so /3. The laser polarization 
n po i can take any direction. We will consider n po i = n 2 , 
the direction perpendicular to the surface, i.e. along the 
symmetry axis of Na§, and one direction orthogonal to it 
with n po i = n T (for the geometry, see section [23t . A laser 
with polarization n z is, of course, an idealization because 
it would correspond to a plane wave running parallel to the 
surface, but it serves to model rather flat impact as com- 
pared to the perpendicular impact of x -polarization. 

From the energy functional, once established, one de- 
rives the static and dynamical equations variationally in a 
standard manner |[T8l . 

2.3 Solution scheme The numerical solution of the 
coupled quantum-classical system proceeds as described 
in lfT9l l20l [39l . The electronic wave functions and the 
spatial fields are represented on a Cartesian grid in three- 
dimensional coordinate space. The numerical box em- 
ployed here has a size of (64 a^f . The spatial derivatives 
are evaluated via fast Fourier transformation. The ground 
state configurations were found by accelerated gradient 
iterations for the electronic wave functions |40J and simu- 
lated annealing for the ions in the cluster and the substrate. 
Propagation is done by the time-splitting method for the 
electronic wave functions ATI and by the velocity Verlet 
algorithm for the classical coordinates of Na + ions and 
MgO or Ar constituents. 

2.4 Gathering angular distributions Electrons 
emitted from the cluster will eventually arrive at the bound- 
aries of the box. In order to suppress re-feed of these elec- 
trons back into the box, we employ absorbing boundary 
conditions 1130114-21 . This is achieved by the mask function 
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for |r| < i? cut 

for Rem < |r| < i? b 

OX 

for Irl > i? box 



(2a) 

where R cut is the cut-off radius outside which absorption 
starts and i?box is the minimal radial distance from the 
origin to the closest point on the boundaries. The Kohn- 
Sham time step actually performed with the time-splitting 
method [41 1 is thus augmented by an absorbing step as 



<P = V 1V <P® 



<p(t + St) = M{r)(p 



(2b) 



where U TV is the unitary propagation operator. Applying 
the mask function A4 to the orbitals gently removes density 
approaching the box boundary and prevents it from being 
reflected. 

To compute the angular distribution of emitted elec- 
trons, the absorbed density is accumulated for each state 
and each (absorbing) grid point as 



/>oo 

A(r) = dt [l-M(r)] U TV 
Jo 



(3a) 



By definition of A4, the field r(r) is non-vanishing only 
in the spherical absorbing zone. The angular distribution of 
emitted electrons is finally gathered by dividing the absorb- 
ing zone into radial segments A v , and integrating -Tj(r) 
over those segments. The Photoelectron Angular Distribu- 
tion (PAD) then becomes 



dN, 



dQ 
df2 



oc 



\A V {6 



W./ r/Kr) ' (3b) 



^ dQ 



(3c) 



where ||t4 w (0, <^)|| denotes the area of the segment A u 
on the surface of a unit sphere. Eq. Oct provides the total 
cross-section, while Eq. Obi the cross-section for emission 
from a specific quantum state (pi. 

We will sometimes present simpler polar distributions, 
averaging the full distributions dN esc (9, <p)/dfi over <f> and 
dividing by the polar volume element sin 9 (see Figure [TJi. 
This is convenient when the variation in </> is weak. How- 
ever, even if the azimuthal distribution may carry interest- 
ing details, particularly for deposited clusters, we will see 
that the reduced view can bring some useful information. 

It is to be noted that we consider here PAD which are 
integrated over all outgoing electron momenta. The state- 
specific PAD carry automatically some information on the 
electron spectra because the dominant emission strength 
goes into the first (multi-)photon peak above continuum 
threshold. More detailed information can be gathered with 
techniques as proposed in 1221 1431 . They also allow to 
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produce a double differential cross-section providing the 
photo-electron spectra in each angular bin separately. An 
explanation of this technique and the much more complex 
analysis of the emerging data is postponed to a later publi- 
cation. 

2.5 Structure of the test cases As test cases, we 
will consider the free Na 8 cluster and Na 8 deposited on 
MgO(OOl) or Ar(OOl) surfaces. Starting point of the laser- 
induced dynamics is a well relaxed structure of each sys- 
tem. These structures had been discussed extensively in 
|. FigureQ]illustrates the structure of Na 8 on MgO(OOl). 




Figure 1 The structure of Na 8 on MgO(OOl). Na ions are 
indicated by black spheres, O ions by white ones, and Mg 
ions by gray ones. The structure of Na 8 on Ar(OOl) is es- 
sentially similar. The bond distance in the two four-fold 
rings of Na 8 is 6.2 ao and the distance between the two 
rings is 5.8 ao- The equilibrium distance between the lower 
cluster plane and the first surface layer is 5 ao. 

The symmetry axis of Na 8 which coincides for the de- 
posited Na 8 with the axis perpendicular to the surface is 
taken as z-axis. The MgO surface is a cut from the cubic 
crystal structure. On the surface, the O and Mg ions are 
arranged in squares where next neighbors are always the 
other species. 

The Na 8 cluster has a highly symmetric configuration 
out of two rings of four ions each, twisted relative to each 
other by 45° to minimize the Coulomb energy. The Na 8 is 
rather rigid and only very little affected by the surface. Free 
Na 8 is very similar to the deposited cluster shown here, 
with bond lengths differing by less than 3%. The same 
holds for Na 8 on Ar(OOl). The electronic cloud of free Na 8 
exhibits close to spherical shape because N e i = 8 electrons 
corresponds to a strong electronic shell closure for Na clus- 
ters fi4l . That changes for the deposited cluster. The sur- 
face destroys the reflection symmetry which, in turn, mixes 
single-electron states of different z-parities. The presence 
of the surface does also shift the single-electron levels and 
the ionization potentials (IP). The relations are sketched in 
figure The IP is 4.3 eV for free Na 8 , 3.8 eV for Na 8 on 
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Figure 2 The single-electron levels of free Na 8 and of Na 8 
deposited on MgO(OOl) or Ar(OOl). The degeneracy of the 
lp x and lp y levels in free Na 8 is indicated by a heavier 
line. The vacuum threshold is at zero energy. The IPs are 
indicated by vertical lines with open arrows and the trans- 
mission barriers by vertical lines with filles arrows. The 
numbers beneath are the IPs or barriers in units of eV. 



MgO(OOl) and 4.2 eV on Ar(OOl). The transmission bar- 
rier for emission from deposited Na 8 into the substrate lies 
at +5.6 eV for MgO(OOl) and +5.4 eV for Ar(OOl): both 
are substantially larger than the IP. The lp x and lp y lev- 
els of free Na 8 are perfectly degenerate, slightly split from 
lp z due to a very small quadrupole deformation of Na 8 . 
The ionic structure of Na 8 deposited on MgO(OOl) hardly 
changes because the magic electron configuration renders 
that cluster very rigid. There is an overall up-shift due to 
core repulsion from MgO. The symmetry breaking due to 
the surface orientation slightly splits the lp x ,y levels, by 
0.204 eV. Much less shifts are seen for Na 8 on Ar(001) and 
the splitting between lp x and lp y , 0.05 eV, is also much 
smaller. That difference correlates to the much lower inter- 
face energy of the Na-Ar system. Indeed, the bottom layer 
of Na 8 is within 5.1 ao from the MgO(OOl) surface, much 
closer to the substrate than in the case of Ar(001) (6.4 ao). 

The optical absorption spectrum of our test cases is 
shown in figure[3]One recognizes the pronounced Mie sur- 
face plasmon resonance ITII3I18II . For free Na 8 , there is one 
clean peak at 2.45 eV almost the same for each mode. The 
surfaces lead to some spectral fragmentation, about 1 .4 e V 
for MgO(OOl) and 0.2 eV for Ar(001). The highly polariz- 
able MgO(OOl) also induces a small down-shift of the reso- 
nance center to 2.3 1 eV while core repulsion outweighs po- 
larization for Ar(001) leading to a small up-shift to 2.54 eV. 
The MgO(OOl) surface leads to a strong fragmentation of 
the z mode essentially due to symmetry breaking l20l . 

In the following studies, we will vary laser frequency 
a>ias and polarization. We work in all cases with the 
same pulse length T pu isc = 60 fs, which corresponds to 
FWHM= 20 fs. The angular distributions in the high- 
intensity regime are always concentrated on forward 
and backward emission along the laser polarization axis, 
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Figure 3 The optical absorption spectrum (dipole re- 
sponse) of free Na§ and of Na§ deposited on MgO(OOl) 
or Ar(OOl). The spectra from the modes in the three ba- 
sic directions are shown separately as indicated. The z axis 
stands for the direction perpendicular to the surface and 
along the symmetry axis of Na 8 . 



whereas a high sensitivity to all details of the excitation 
is found in the regime of weak perturbations 1231 . The 
present study thus aims at staying close to the perturba- 
tive regime to explore the rich variety of distributions. The 
intensity should be low enough for a perturbative regime 
being valid |22 26 ,45 1, but also sufficiently high to provide 
signals safely above numerical noise. We found an inten- 
sity of / = 10 9 W/cm 2 to be a good compromise and we 
used that value in most of the test cases. Use of a different 
value will be clearly indicated. 

3 Brief survey of free Na 8 Detailed studies of the 
PAD of free clusters will be discussed in a forthcoming 
publication. We briefly summarize here the results. 

Variation of the laser intensity shows almost constant 
pattern of the PAD throughout the regime of moderate in- 
tensities. For further increasing intensities, the structures 
change steadily towards simple forward-backward scatter- 
ing. This is related to the transition from the frequency- 
dominated regime of moderate intensities to the field dom- 
inated highly non-linear regime |27l l30i[46ll47l . 

Variation of frequency leads to dramatic changes in the 
angular distributions. They vary amongst three typical pat- 
terns. We have chosen three frequencies to have an ex- 
ample for each type. Results are shown in figure |4] The 
most common case is forward scattering (along the axis 
of polarization, 9 = 0,tt), seen here in the lowest panel 
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Figure 4 Gray scale plot of angular distributions for free 
Na§ excited with three different laser frequencies, as indi- 
cated. The gray scale is used in arbitrary units. High emis- 
sion is signified by white and no emission correlates to 
deep black. The frequencies had been selected to display 
the three different patterns which could be found. The laser 
is polarized along the z-axis, the symmetry axis of Na 8 . 

for wi as = 2.5 eV. Sometimes one encounters "diagonal 
scattering" where electron emission is concentrated on a 
double cone around angle 7r/4 (with respect of the z-axis), 
see upper panel for cji as = 5.4 eV. Sideward scattering as 
exemplified in the middle panel, is observed occasionally, 
here for wi as = 3.0 eV. Note that the spectral relations of 
these three cases are very different. For 5.4 eV, we have 
a one-photon process for the lp states but a (much sup- 
pressed) two-photon process from the Is state. For 3 eV, 
we have a two-photon process for both, the lp and the Is 
shell. And for 2.5 eV, we have a two-photon emission from 
the lp while the Is shell, again, requires one more photon. 

The strong frequency dependence of the PADs can be 
related to the fact that the strongest emitting single-electron 
state depends sensitively on the relationship between laser 
frequency and IP, and that the emission from each state 
looks very differently. This is demonstrated by analyzing 
the PADs of individual single-electron states. Experimen- 



Copyright line will be provided by the publisher 



6 



Matthias Bar et al.: Angular distributions of electrons from Nag 



tally, this is achieved by measuring the PAD simultane- 
ously with the photoelectron spectrum (PES) I8l [l0ll48l . 
The energy selection by PES allows one to associate the 
PAD to the single-electron states from which these were 
produced. Theoretical calculations have the separate in- 
formation trivially at hand, as seen in Eq. Oct . Figure [5] 




90° 180 270° 360° 



4> 

Figure 5 State selective angular distribution 
dN OSCt i(9,(f>)/df2 of electrons which were emitted 
from single-electron state tpi of free Nag as indicated, 
given in arbitrary units. The laser was polarized along the 
z axis and had a frequency of cji as = 2.5 eV. 

shows the state-specific PAD of free Nag for the (hypothet- 
ical) case that the clusters symmetry axis is aligned with 
the laser polarization. The separate distributions are much 
more structured and show their maxima at different an- 
gles. Particularly noteworthy are the pattern from the lp x 
and lp y states. Both have their emission maxima at po- 
lar angle 9 = 45° and 135° and both show pronounced 
structures in azimuthal angle (f>; the azimuthal pattern is 
shifted by 90° between lp x and lp y . That feature will 
play a role again for deposited Nag, see section [4] The 
maxima at 6 = 45° and 135° can be qualitatively under- 
stood in a picture using wave functions of a spherical mean 
field. The lp xy have then the angular wave function in 



terms of the spherical harmonics Y±. The dipole excita- 
tion comes with angular distribution Y®. For lp x , the emit- 
ted wave is driven by the product + Y{~ )Y]°| 2 oc 
sin 2 (28) sin 2 </>, having maxima at <f> = 0° and 90°, and for 
lp y by - Yf^Yfl 2 oc sin 2 (26») cos 2 with max- 
ima at <j> = 45° and 135°, while both waves are maximized 
at 9 = 45° and 135°. This picture agrees with the obtained 
patterns in figure|5] Note that the total cross-section of free 
Na§ does not exhibit any structure in <fi as the lp x and lp y 
distributions add up to a constant value in <f> (besides a 
faint perturbation by the ions). A similar finding holds for 
the energy-resolved distributions as the states lp x and lp y 
have exactly the same single-electron energy and thus are 
added up with equal weight. We also separately explored 
the frequency dependence of the distributions for emission 
from each single-electron state. They turned out to depend 
much less on the frequency than the total PAD. We refrain 
from presenting these results in detail as this would require 
quite a bit space. 
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Figure 6 Angular distributions, averaged over orienta- 
tions, along polar angle 8 for emission from the different 
single-electron states of free Na 8 . The laser frequency is 

w las = 5.44 eV. 

Considering free clusters with known orientation al- 
lows one to investigate characteristic structural features 
and to provide basic information for deposited clusters to 
be discussed later on. An alignment of clusters in the gas 
phase remains an open experimental problem ("orientation 
burning" may be an option ||49l ). Present-days experiments 
in the gas phase measure an ensemble of clusters where all 
orientations are equally distributed. A simulation of that 
situation requires one to calculate an ensemble of orien- 
tations and to average the emerging angular distributions. 
The procedure is straightforward. We computed angular 
distributions on a grid of orientations with grid spacing of 
22.5° in Euler angles 9 and <f>, ignoring ip due to the near 
axial symmetry of Nag, and added them with appropriate 
geometrical weights. The orientation averaging wipes out 
the sub-structures along <f> direction. Figure|6]shows the re- 
sulting (orientation averaged) distributions along polar an- 
gle 9. Note that we have also averaged over lp x and lp y 
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Figure 7 Angular distributions for emission from Nag 
on MgO(OOl) for three different laser frequencies as indi- 
cated. The laser is polarized along the z-axis, perpendicular 
to the surface. 



distributions since these two states are energetically degen- 
erate and could not be discriminated by PES. The orienta- 
tion averaging yields much different and simpler structures 
than those seen in the case of known orientation, see fig- 
ure [5] The dominance of purely forward-backward struc- 
ture holds only for the small system Nag. Larger clusters 
show richer structures 



4 Emission patterns for deposited Na 8 
4.1 Deposition on MgO(001) As pointed out in sec- 
tion |3] the detailed analysis of free PAD leads to remark- 
able insights into the structures of clusters. However, the 
orientation problem inherent to clusters in the gas phase 
blurs some of the detailed information. It is thus interesting 
to consider the complementing case of deposited clusters 
where the orientation is well defined. On the other hand, 
the presence of the substrate complicates the PAD. Analy- 
sis of this effect is thus a key issue. 

Figure [7] displays angular distributions of Nag on 
MgO(OOl) for three different frequencies. Again, the fre- 
quencies were selected to show three different and typical 
emission patterns. The sequence is comparable to the case 
of free Nag discussed above, but having somewhat lower 
values to accommodate for the lower IP (see figure [2]i and 
stay below the transmission threshold. The substrate ob- 



viously has a very strong influence such that the patterns 
are quite different from those of free Nag. The dominant 
forward emission (8 = 0°) is still observed; backward 
emission (6 = 180°) is, of course, totally suppressed by 
the presence of the insulating substrate. More surprising is 
the appearance of a strong azimuthal dependence. Recall 
the pronounced azimuthal structures for emission from 
the lp x . y levels of free (aligned) clusters, see figure [5] In 
the total or energy-resolved cross-sections of free clusters, 
these structures are wiped out by the degeneracy, as each of 
the two levels contributes equally to the emission strength, 
see the discussion of figure [5] This lp x -\p y degeneracy 
is now split by the surface leading to an energy difference 
between lp x and lp y state of 0.2 eV. Because these states 
are close to the emission threshold, the small energy differ- 
ence has a large effect on the relative emission strengths. 
One of the two states dominates emission and its pattern 
affects the total distribution. Comparison to figure[5]makes 
it clear that the dominant state is lp x for oji as = 2.6 eV 
(figure |7] lower panel) and lp y for wi as = 3.7 eV (figure 
[7] middle panel). There is also some effect on the po- 
lar distribution. The emission cones are neither perfectly 
aligned along 6 = 90°, as in a sidewards emitting case, nor 
along 6 = 45°, as was the case for lp x . y orbitals of free 
Nag. The composition of emission strengths from the four 
single electron levels is different for the deposited case 
as compared to free Nag. Furthermore, the polarization 
attraction bends the outwards cones towards the substrate. 
Both effects together widen the emission cone such that 
the angle of inclination of the emitted electrons relative to 
the surface decreases. 

This example demonstrates that the detailed structure 
of the angular distribution, in the polar as well as in the az- 
imuthal angles, sensitively depends on the laser frequency. 
There is much more structure compared to the emission 
spectra of free clusters, thus more information, worth to be 
unraveled. That information, however, is masked by the in- 
volved interplay of cluster and surface which inhibits sim- 
ple one-to-one correspondences and requires a careful cor- 
relation analysis. 

4.2 Comparison with argon substrate We also in- 
vestigated the case for another combination, namely Nag 
on Ar(001). Frozen Ar is an insulator as MgO is. We thus 
expected, as in the case of Nag on MgO(OOl), heavily 
modified angular distributions compared with the free Nag 
case. Ar is, however, a Van-der-Waals bound system with a 
smaller surface corrugation and smaller initial polarization 
fields, in contrast to the ionic crystal MgO. Figure|8]shows 
angular distributions for Na 8 on Ar(001) at two frequen- 
cies of the laser field. The effects are very similar to the 
case of MgO. Backward scattering is totally suppressed, 
forward scattering is accordingly enhanced, there appear 
pronounced azimuthal structures, and the distributions de- 
pend sensitively on the frequency. With 0.14 eV, the split- 
ting of the lp x . y levels due to Ar(001) surface is similar to 
that at MgO, which is probably responsible for the similar- 
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Figure 8 Angular distributions for emission from Na 8 on 
Ar(OOl) for two laser frequencies as indicated. The laser is 
polarized along the z-axis, perpendicular to the surface. 

ities of the emission patterns between Ar and MgO. There 
is, however, one difference to the case of MgO. The side- 
ward cone (lower panel of figure [8) is not so close to the 
surface and points towards a polar angle of 6 = 45° as in 
the case of free Nag. That is probably due to the smaller 
surface polarization for Ar. 

4.3 Frequency and intensity dependences In this 
section we explore more systematically the influence of 
laser frequency and intensity on the angular distributions 
for deposited Nag. To simplify the graphical represen- 
tation, we consider </>-averaged distribution along polar 
angle 8. Figure|9]shows azimuthally averaged distributions 
for emission from Nag on MgO(OOl). The upper panel 
presents the variation with u\ as for moderate intensity. 
The prevailing effect is the strong suppression of back- 
ward emission through the substrate. In spite of reflection- 
enhanced surface dominance, we see again the variation of 
pattern with frequency qualitatively similar to the case of 
free Nag (see figure [4j. The pattern is, however, more in- 
volved than what would emerge from a simple rescattering 
of the reflected electrons to angle 6 — > 180° — 6. The sur- 
face attraction shifts the former diagonal maximum from 
45° towards the substrate and produces sizeable emission 
under rather flat angles 9 w 90°. A first glimpse of trans- 
mission into the substrate can be seen for wi as = 5.44 eV. 
That is very close to the transmission threshold for Na 8 on 
Mg(001) at about 5.58 eV, see figure0 

The lower panel of figure [9] shows the changes with 
increasing intensity for fixed frequency. We see again the 
expected typical increasing focus towards forward emis- 
sion. A somewhat surprising and most interesting effect is 
that backward emission gathers rather large strength with 
further increasing intensity. The high field strengths make 




polar angle e [degrees] 

Figure 9 Top : Azimuthally averaged angular distribu- 
tion (arbitrary units) of electrons emitted from Nag on 
MgO(OOl) irradiated by laser pulses with intensity / = 10 9 
W/cm 2 and varying frequency. The laser is polarized along 
the z-axis, perpendicular to the surface. Bottom : Distribu- 
tion for fixed frequency wi as = 4.76 eV with varying laser 
intensity. 



two- and more-photon processes competitive which, in turn 
overrules the frequency counting. However, one has to be 
cautious with interpreting these processes where substan- 
tial amounts of electron charge are penetrating into the sub- 
strate. The present QM/MM model is not set up for that 
situation and likely is no longer accurate at a quantitative 
level. 

Figure [10] shows some azimuthally averaged distribu- 
tions for Na 8 on Ar(001). The general trends are very sim- 
ilar to the case of MgO, see for example the upper panel of 
figure[9] We see again the strong forward dominance due to 
reflection from the substrate and the strong dependence of 
the distribution on laser frequency as was observed before 
for all cases. There is a minor difference in that the chances 
for transmission into the substrate are even slightly lower. 

5 Conclusion We have explored from a theoretical 
perspective angular distributions of electrons emitted from 
Na 8 clusters, free and deposited on Ar(001) or MgO(OOl) 
surfaces. Thereby we concentrated on direct electron emis- 
sion which takes place within few fs after laser excitation 
and which is the dominant process for short laser pulses (up 
to about 50 fs). Several laser frequencies were used, higher 
ones above the threshold for one-photon emission, lower 
frequencies in the two-photon regime close to the Mie 
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Figure 10 Azimuthally averaged angular distributions (ar- 
bitrary units) for Nag on Ar(001) for intensity I = 10 9 
W/cm 2 and three different laser frequencies, as indicated. 
The laser is polarized along the z-axis, perpendicular to the 
surface. 

plasmon resonance, and cases safely off-resonance. The 
laser intensity was in most cases / = 10 9 W/cm 2 , large 
enough to overcome numerical noise, but safely below the 
regime of violent and highly non-linear processes. The 
numerical simulations employed a hierarchical quantum- 
mechanical-molecular-mechanical (QM/MM) model ap- 
proach. The cluster electrons were described by time-de- 
pendent density-functional theory (TDDFT) at the level of 
the local-density approximation (LDA). A self-interaction 
correction has been added to put the single-electron levels 
to the appropriate energy relative to the particle continuum. 
The cluster ions were propagated by classical molecular 
dynamics. Similarly, the Ar(001) or MgO(OOl) substrates 
were treated as classical systems with atomic positions 
and dipole polarizability as dynamical degrees of freedom. 
The TDDFT equations were solved on a three-dimensional 
coordinate-space grid without any symmetry restriction. 
Absorbing boundaries were applied and the angular distri- 
butions were obtained by recording the electron absorption 
at each grid point in the absorbing zone. 

We first briefly explored free Na 8 . The state-specific 
distributions for Na 8 aligned with laser polarization show 
pronounced patterns in both angles, 9 and <f>, which ex- 
hibit clear footprints of the emitting state. Simulating the 
experimental situation which deals with orientation aver- 
aged ensembles renders the distributions independent of 
the azimuthal angle (f> and reduces structures in the polar 
angles 9. A study of averaging effects, size-, frequency-, 
and intensity-dependence for free clusters will follow in an 
upcoming publication. 

The present study focused on angular distributions 
from Na 8 deposited on MgO(OOl) and Ar(001). The at- 
tachment to the surface provides a well defined cluster 
orientation and it adds substantial perturbations from the 
surface interaction. The large band gaps of both materials 
manifested in high transmission barriers, almost totally 
suppress backward emission, focusing electrons into a 



forward cone. On the other hand, the long-range polariza- 
tion attraction can bend sideward flowing electrons down 
towards the surface, thus widening the opening angle of 
the emission cone. Polarization is strong for MgO(OOl) 
which induces a sizeable trend to emission parallel to the 
surface. Repulsion dominates for Ar(001) which stabilizes 
the forward cone. The surfaces do also break the four- 
fold symmetry (C4) of Nag which, in turn, removes the 
degeneracy of the two single-electron levels with nodes 
orthogonal to the symmetry axis, and thus inhibits the 
compensation of <\> dependence which was observed for 
(oriented) free Na 8 . Pronounced azimuthal structures are 
found for the distributions of the deposited cases. 

Altogether, the computational study revealed interest- 
ing structures which are, however, somewhat complicated 
by surface (particularly the electronic and geometric struc- 
ture of the clusters) contained in the calculated photoelec- 
tron angular distributions and their energy dependence. 

Acknowledgements This work was supported by the 
Deutsche Forschungsgemeinschaft (RE 322/10, RO 293/27), 
Fonds der Chemischen Industrie (Germany), a Bessel-Humboldt 
prize, and a Gay-Lussac prize. 

References 

[1] M. Brack, Rev. Mod. Phys 65, 677 (1993). 

[2] W. A. de Heer, Rev. Mod. Phys. 65, 611 (1993). 

[3] U. Kreibig and M. Vollmer, Optical Properties of Metal 

Clusters (Springer Series in Materials Science, vol 25, 

1993). 

[4] H. Haberland (ed.), Clusters of Atoms and Molecules 1- 
Theory, Experiment, and Clusters of Atoms (Springer Se- 
ries in Chemical Physics, vol 52, Berlin, 1994). 

[5] K. M. McHugh, J. G. Eaton, G. H. Lee, H. W. Sarkas, L. H. 
Kidder, J. T. Snodgrass, M. R. Manaa, and K. H. Bowen, J. 
Chem. Phys. 91, 3792 (1989). 

[6] D. L. Lichtenberger, K. W. Nebesny, C. D. Ray, D. R. Huff- 
man, and L. D. Lamb, Chem. Phys. Lett. 176, 203 (1991). 

[7] J. Pinare, B. Baguenard, C. Bordas, and M. Broyer, Eur. 
Phys. J. D 9, 21-24(1999). 

[8] B. Baguenard, J. C. Pinar, C. Bordas, and M. Broyer, Phys. 
Rev. A 63, 023204(2001). 

[9] J. R. R. Verlet, A. E. Bragg, A. Kammrath, O. Chesh- 
novsky, and D.M. Neumark, J. Chem. Phys. 121, 10015- 
10025 (2004). 

[10] O. Kostko, C. Bartels, J. Schwobel, C. Hock, and B. von 
Issendorff, J. Phys. : Conf. Ser. 88, 012034 (2007). 

[11] S. Skruszewicz, J. Passig, A. Przystawik, 
J. Tiggesbaumker, and K. H. Meiwes-Broer, to be 
published (2009). 

[12] C. Bartels, C. Hok, J. Huwer, R. Kuhnen, J. Schwobel, and 
B. von Issendorff, Science 323, 132 (2009). 

[13] K. H. Meiwes-Broer (ed.), Metal clusters at surfaces 
(Springer, Berlin, 2000). 

[14] K. H. Meiwes-Broer (ed.), Clusters at Surfaces: Electronic 
Properties and Magnetism (Applied Phys. A 82, special is- 
sue, 2006). 

[15] M. Moseler, H. Hakkinen, and U. Landman, Phys. Rev. 
Lett. 89, 176103 (2002). 



Copyright line will be provided by the publisher 



10 



Matthias Bar et al.: Angular distributions of electrons from Nag 



[16] V. V. Semenikhina, A. G. Lyalin, A. V. Solov'yov, and 
W. Greiner, JETP 106, 678 (2008). 

[17] R. M. Dreizler and E. K. U. Gross, Density Functional The- 
ory: An Approach to the Quantum Many-Body Problem 
(Springer- Verlag, Berlin, 1990). 

[18] P. G. Reinhard and E. Suraud, Introduction to Cluster Dy- 
namics (Wiley, New York, 2003). 

[19] F. Fehrer, P. G. Reinhard, E. Suraud, E. Giglio, B. Gervais, 
and A. Ipatov, Appl. Phys. A 82, 151 (2005). 

[20] M. Bar, L. V. Moskaleva, M. Winkler, P. G. Reinhard, 
N. Rosch, and E. Suraud, Eur. Phys. J. D 45, 507 (2007). 

[21] E. Giglio, P. G. Reinhard, and E. Suraud, Phys. Rev. A 67, 
43202 (2003). 

[22] A. Pohl, PhD thesis, Friedrich-Alexander-Universitat, Er- 

langen/Niirnberg, 2003. 
[23] A. Pohl, P. G. Reinhard, and E. Suraud, Phys. Rev. A 70, 

023202 (2004). 

[24]E.E.B. Campbell, K. Hansen, K. Hoffmann, G. Korn, 

M. Tchaplyguine, M. Wittmann, and I. V. Hertel, Phys. 

Rev. Lett. 84, 2128 (2000). 
[25] R. Schlipper, R. Kusche, B. von Issendorff, and H. Haber- 

land, Appl. Phys. A 72, 255 (2001). 
[26] A. Pohl, P. G. Reinhard, and E. Suraud, J. Phys. B 37, 3301 

(2004). 

[27] T. Fennel, K. H. Meiwes-Broer, J. Tiggesbaumker, P. M. 
Dinh, P. G. Reinhard, and E. Suraud, Rev. Mod. Phys., in 
press (2010). 

[28] M. Biir, PhD thesis, Friedrich-Alexander-Universitat, Er- 

langen/Niirnberg, 2008. 
[29] P.M. Dinh, P.G. Reinhard, and E. Suraud, Phys. Rep., in 

press (2009), http://arxiv.org/abs/0903.1004vl. 
[30] F. Calvayrac, P. G. Reinhard, E. Suraud, and C. A. Ullrich, 

Phys. Rep. 337, 493 (2000). 
[31] V. A. Nasluzov, K. Neyman, U. Birkenheuer, and N. Rosch, 

J. Chem. Phys. 115, 17 (2001). 
[32] F. Fehrer, M. Mundt, P. G. Reinhard, and E. Suraud, Ann. 

Phys. (Leipzig) 14, 411 (2005). 
[33] P. M. Dinh, F. Fehrer, P. G. Reinhard, and E. Suraud, Eur. 

Phys. J. D 45, 415(2007). 
[34] P. M. Dinh, F. Fehrer, P. G. Reinhard, and E. Suraud, Surf. 

Sci. 602, 2699 (2008). 
[35] C. Legrand, E. Suraud, and P. G. Reinhard, J. Phys. B 35, 

1115 (2002). 

[36] G. R. Ahmadi, J. Almlof, and J. Roegen, Chem. Phys. 199, 
33 (1995). 

[37] F. Duplaa and F. Spiegelmann, J. Chem. Phys. 105, 1492 
(1996). 

[38] M. Gross and F. Spiegelmann, J. Chem. Phys. 108, 4148 
(1998). 

[39] F. Fehrer, P. G. Reinhard, and E. Suraud, Appl. Phys. A 82, 
145 (2005). 

[40] V. Blum, G. Lauritsch, J. A. Maruhn, and P. G. Reinhard, J. 

Comp. Phys 100, 364 (1992). 
[41] M. D. Feit, J. A. Fleck, and A. Steiger, J. Comp. Phys. 47, 

412 (1982). 

[42] C. A. Ullrich, J. Mol. Struct. (THEOCHEM) 501-502, 315 

(2000). 

[43] A. Pohl, P. G. Reinhard, and E. Suraud, Phys. Rev. Lett. 84, 
5090 (2000). 



[44] S. Bjornholm and J. Borggreen, Phil. Mag. 79, 1321 
(1999). 

[45] F. H. M. Faisal, Theory of Multiphoton Processes (Plenum 

Press, New York, 1987). 
[46] P.G. Reinhard, F. Calvayrac, C. Kohl, S. Kummel, 

E. Suraud, C. A. Ullrich, and M. Brack, Eur. Phys. J. D 

9, 111 (1999). 

[47] G. Zwicknagel, C. Toepffer, and P. G. Reinhard, Phys. Rep. 

309, 117 (1999). 
[48] O. Kostko, B. Huber, M. Moseler, and B. von Issendorff, 

Phys. Rev. Lett. 98, 043401 (2007). 
[49] P. G. Reinhard and E. Suraud, Eur. Phys. J. D 34, 145 

(2005). 

[50] M. Baer, P. M. Dinh, L. V. Moskaleva, P. G. Reinhard, 
N. Rosch, and E. Suraud (2009), in preparation. 

[51] C. Barrels, PhD thesis, Albert-Ludwigs-Universitat, 
Freiburg im Breisgau, 2008. 



Copyright line will be provided by the publisher 



